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Introduction 



The s-d exchange model, often referred to as the Kondo model, has been extensively discussed since 
it was first used by Kondo [1] to calculate the effects of impurities on the resistivity of a metal. Various 
approaches have been applied to study this model, including renormalization group methods [2], the Bcthe 
ansatz [3]-[9] and conformal field theory [10]-[12] . Traditionally the Kondo model contains a single coupling 
J. Allowing the exchange couplings to differ gives the anisotropic Kondo model. Regardless of whether there 
is an anisotropy or not, the Kondo model is a massless integrable system for arbitrary spin [13], [14]. 

In this paper we present two massive generalizations of the anisotropic spin 1/2 Kondo model and 
discuss their integrability. These models are bosonic boundary field theories related to the original fermionic 
system via bosonization. Such massive theories have possible applications to ID impurity problems with an 
excitation gap. An example is a magnetic impurity in a superconductor, where the mass represents the BCS 
energy gap [15]. 

The motivation for considering our particular models came as follows. The massive boundary sine- 
Gordon model (MBSG) is a boundary field theory of a single boson <f) with a bulk interaction in the action 
G J dxdt cos((3(j)) and a boundary interaction A / dtcos(f3(4> — 4> )/2). The model is known to be integrable 
for arbitrary G, A and 4> . The massless limit G — > is well defined, which defines the massless boundary 
sine-Gordon model [16]. On the other hand, the usual anisotropic Kondo model after bosonization has 
a boundary sine-Gordon like interaction but with the inclusion of spin operators S± at the boundary: 
Xjdt S+e 2 ^/ 2 + S'-e~ 1 ' 3 ^/ 2 , and the bulk is just a free massless scalar field. If one compares the massless 
boundary sine-Gordon with the Kondo model, one finds that the same sine-Gordon spectrum of particles 
diagonalizes the boundary interaction; the difference is in the detailed form of the scattering matrices for 
reflection off the boundary [17]. So the question naturally arises whether it is possible to add a bulk sine- 
Gordon term to the Kondo theory and preserve the integrability, in analogy with the massive boundary 
sine-Gordon theory. 

If the theories we define in this paper are indeed integrable this implies that the usual bulk sine-Gordon 
integrals of motion are not spoiled by the boundary interactions. In this work we do not establish the 
integrability of our models by studying these integrals of motion. Rather we take the following approach to 
studying the problem. In [18] the most general solution to the boundary Yang-Baxter equation corresponding 
to a massive sine-Gordon spectrum and no boundary degrees of freedom was described and related to the 
massive boundary sine-Gordon model. The known scattering description of the Kondo model also does 
not require boundary degrees of freedom for spin 1/2 due to screening. In the massive theory one might 
expect a competition between the gap effects and the screening, since the screened state is a property of a 
new infra-red fixed point, but the existence of the gap introduces an infra-red cutoff, and this may imply 
that this infra-red fixed point is never reached. However non-trivial renormalization group beta functions 
only occur at the isotropic point. Our analysis of the free fermion point indeed indicates that a variant of 
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screening does occur in the massive case. We thus assume this continues to hold for the massive model, 
namely that the scattering description does not require any boundary degrees of freedom, at least in some 
sector. Then if the massive version of the Kondo model is integrable, since it involves a bulk sine-Gordon 
theory, its reflection S-matrix off the boundary must be contained in the general solution found by Ghoshal 
and Zamolodchikov (GZ) [18] up to overall scalar (CDD) factors. If one simply adds a bulk sine-Gordon 
term to the usual Kondo model, and specializes to the free fermion point, then the reflection amplitudes can 
be computed explicitly and it is found that they do not coincide with those of GZ, even up to CDD factors. 
However we found that a slight modification of the boundary interaction does match onto the GZ solution 
and it is this model for which we conjecture an extension away from the free fermion point. This latter model 
requires that the boundary and bulk couplings not be independent. The first model on the other hand does 
not have the constraints of the boundary Yang-Baxter equation at the so-called reflectionless points, so in 
this situation we also conjecture an exact reflection S-matrix. 

In the massless cases the main distinction between the Kondo and boundary sine-Gordon theories is that 
sine-Gordon has a flow between free boundary conditions in the ultra-violet and fixed boundary conditions in 
the infra-red, whereas the Kondo model maintains a free boundary condition throughout. Our conjectured 
reflection S-matrix for the modified massive Kondo model is in accordance with this, in that it corresponds 
to the reflection S-matrix for the boundary sine-Gordon theory at the free boundary condition times a CDD 
factor. 

We present our results as follows. In section 1 we define the two massive versions of the Kondo model 
that we consider. The first, simply referred to as massive Kondo (MK), is just a sine-Gordon model with 
the usual Kondo interaction at the boundary. The second, the modified massive Kondo model (mMK), has 
a slightly different boundary interaction. More importantly, in the mMK model the boundary coupling A is 
not independent of the bulk coupling G, but satisfies G a A 2 . This is reminiscent of what happens in the 
boundary Toda theories [19]. In section 2 we study both models at the free fermion point, and relate them 
both to two decoupled Ising models in appropriate boundary magnetic fields. In section 3 we compare with 
the GZ solution and propose reflection S-matrices away from the free fermion point. Section 4 contains our 
results on boundary bound states for the massive models. An appendix reviews the mapping of the Kondo 
model to a bosonic boundary field theory. 

1. The Field Theory Models 

The anisotropic Kondo Hamiltonian is 

hK = T, + y Sz E 4>^<w + # E (^4>*)- c ^ + ^4> a )- c ^) . 

her kk'aa' kk'aa' 

(l.i) 

where ct^ are conduction electron creation operators with wave vector k and spin a (up \ or down j), a x ' v ' z 
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are the Pauli matrices, s is the impurity spin operator (spin 1/2) at r = and J z and J± are the exchange 
couplings. In the isotropic Kondo model J z = J±. For s-wave scattering near the fermi surface, the spin 
sector of (1.1) can be mapped onto the following bosonic boundary field theory 

H K = \f dx ((d t ct>) 2 + (d x 4>) 2 ) + A (s+e^' 2 + S_ e -WW/ 2 ) , (i.2) 



where A is the coupling and (3 determines the anisotropy. The isotropic point corresponds to (3 = \/8ir. The 
series of steps leading to (1.2) are outlined in the appendix. 

The boundary field theory (1.2) is our starting point, which we will simply refer to as the Kondo 
model. The parameters A and (3 are considered to be arbitrary parameters, independent of their relation 
to the couplings J z and J± (see appendix). This model is known to be integrable for any spin. However, 
integrability requires that the matrices Si form a spin j representation of the g-deformed quantum algebra 
su(2) q [13], where q — exp(i/3 2 /8). The su(2) q relations are 

[S Z ,S ± ]=±2S ± , [S + ,S.]= qSZ q ~^ Z . (1.3) 

For the isotropic case, q = — 1 and (1.3) reduces to the usual su(2) algebra. For spin 1/2, the difference 
between su(2) q and su(2) is not important since the Pauli matrices also form a representation of su{2) q for 
arbitrary q. 

1.1. The Massive Kondo Models 

We are interested in a massive generalization of (1.2), which in the massless limit reproduces the Kondo 
behaviour. In particular, the massless limit should give the spin 1/2 Kondo scattering matrix. This has led 
us to consider two models. 

The first model, the MK model, has the Hamiltonian 

H MK = \f dx ((d^) 2 + (d x 4>? - Gcos(/ty)) + A (s+e^/ 2 + S_ e -W(°>/ 2 ) . (1.4) 

The bulk part is simply the sine-Gordon model. This model was studied in [20] at the free fermion point 
((3 = VItt) and the reflection S-matrix was calculated. We will re-derive the free fermion reflection amplitudes 
and discuss their extension to arbitrary (3. It is known that the bulk sine-Gordon spectrum also diagonalizes 
the boundary interaction [18], thus (1.4) is perhaps the most obvious massive generalization. The bulk 
spectrum consists of solitons for < (3 < V&k, and for (3 < \I~A-k there are also breathers. (For (3 > V8tt 
the bulk energy is unbounded.) The soliton and breather masses are functions of G and (3. We will argue that 
unlike the Kondo model, (1.4) is not integrable for arbitrary /?, but only at certain so-called reflectionlcss 
points. 

This brings us to the second model, the mMK model, with the Hamiltonian 

H mMK = ^ dx + {dx4)) 2_ Gcosm ^ 

-iX (cos |(^(0) - 4>)S+ - S- cos |(^(0) - 4>*)\ , (1.5a) 



where 

0--|(f-^o), (1.5b) 
and the coupling A is related to the bulk coupling G 

A 2 cx G, A = V2M at /? = V^r. (1.5c) 

Here 0o is an arbitrary real parameter and M is the fermion (soliton) mass. Let us compare this with the 
MBSG model 

H MBSG = l -f dx {(d t <pf + {8 X 4>) 2 -Gco S (f3^))+\ C os^{4>(Q)-M- (1-6) 

^ J — oo ^ 

The MBSG model is integrable for all values of A, (j) and G, and all can be considered as independent [18]. 
In the mMK model there is only one free coupling (j) , with A being a function of the bulk coupling G, i.e. 
not a free parameter. We believe that the mMK model is integrable for all values of (3 and conjecture a 
soliton reflection S-matrix in section 3, which in the massless limit recovers the Kondo model. Note that for 
both massive models, the scaling dimension of the bulk coupling G is twice that of the boundary coupling A. 
Because of the cosines and exponentials, these couplings have anomalous dimensions, resulting in a scaling 
dimension of 2 - (3 2 /Air for G and 1 - (3 2 /8n for A. 



2. Boundary Reflection S-Matrices at the Free Fermion Point 

In this section we derive the reflection S-matrices at the free fermion point starting from the action. 
A similar procedure was used in [21] for the MBSG model. For the MK model this was first done in [20]. 
Here we repeat that calculation in a more general context, and by making a change of basis give a clearer 
understanding into the structure of the reflection S-matrix. 

The bulk action for both models is 



(2.1) 



-i o OO p u 

Sbuik = - / dt dx ({d t <t>) 2 - (d x <f>) 2 + Gcos(/ty)) . 

^ J — oo J — oo 

In the massless limit, the scalar field can be separated into its right and left moving components 

<p(z + ,z-) = <p{z+)+7p{z-), z ± =t±x, (2.2) 
in terms of which the Neumann boundary condition d x (f) = at x = becomes 

d z+ v(z + )=d z -lp( Z -), (x = 0). (2.3) 

This implies 

p(s+)= (* = 0), (2.4) 

\/47T 
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where a is a constant related to the zero modes of 4> [21], [22]. The massless boson bulk theory is equivalent 
to a massless Dirac theory through bosonization/fermionization. Let ip±(z + ) and ip±(z~) respectively be the 
left and right chiral components of the Dirac fermion with charge ±1. Then the bosonization relations 
read 

ip±(z+) = e ±*v^P(* + ) ; $ ± { z -) = e W4^(0 ; ^_ = ^ + )t ; ^_ = (^ + )t_ (2.5) 
In terms of the fermions the (free) boundary condition becomes 

V>± - e Tl > T , (x = 0) (2.6) 

which breaks the charge symmetry. 

For the massive case, the bulk theory (2.1) is equivalent to a free massive Dirac theory at /3 = 
[23], [24]. The fields (tp± 7 ip±) are now no longer chiral. These can still be related to <f) as in (2.5), but now ip 
and Tp, with <p = f + ty, arc not simply the left and right moving components [25] . The fermions ip± and tp± 
correspond to the solitons in the sine-Gordon model. Considering the massive system as a perturbation of 
the massless one, we use the same boundary conditions (2.4) and (2.6). The massive fermionic action which 
enforces (2.6) takes the form 

S = S'buik + Sbc (2-7) 

Sbulk — Skin + S mass (2.8a) 

/oo [>0 
dt / dx [ii>-(d t - d x )xjj+ + ixjj+{d t - 9 x )v_ + i$-(d t + d x )^ + + hp + {d t + d x y$_] (2.8b) 
-00 J — 00 

/oo 

dt dxM(^J^ + ^_^ + ) (2.8c) 
-OO J —00 

/OO 
dt (e 4CT V+?+ + e~ ifT i>-^_) , (2.9) 
-00 

where M is the soliton mass. The action S'buik is simply a Dirac action. S\, c is required to produce the 
boundary conditions. In varying S with respect to all the fields, the boundary terms give (2.6). Note that 
the i/>'s are different from the original fermions occurring in the Kondo model. We went from (1.1) to the 
bosonic boundary field theory for cf) (1.2) by bosonizing (see appendix), and now we have (i) added a mass 
term and (ii) re-fermionized, going from <j> to a new set of Dirac fermions. The mass term (2.8c) is not a 
mass term for the original Kondo fermions. (In the Kondo model a similar transformation gives the resonant 
model [26], [33].) The action (2.7) is common to both the MK and mMK model. 

We now specify the interaction with the spin operators S±. It will be convenient to consider a more 
general boundary interaction, with the MK and mMK interactions being special cases. Consider the following 
action 



-Ax dt (e^' 2 S+ + S_ e -WW/ 2 ) + A 2 r dt (e-^/ 2 S+ + S_ e W( )/ 2 ) , (2. 

J — oo J — oo ^ 
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where Ai and A 2 are couplings which can be different. The MK model is obtained by setting A2 = and 
Ai = A. For the mMK model we need to take X 1 = Ae*> and A 2 = Ae^ . Using (2.4) and (2.5) the 
interaction term (2.10) can be written as 

A f°° 

Sint = --z- / dt[(ip + a-+ip_a+)S+ + S-{a + ip- + a-ip + )] 
z J-00 

A f°° 

+ -1 / dt[(a+il>-+a-il> + )S + + S-(il; + a-+il)_a + )], (2.11) 
2 J-00 

where a± are the zero mode contributions 

a+ = e- ia/2 , a- = (a+)t = e la/2 , (2.12) 

and we have taken the average in fermionizing 

e +iV**(0) = l(ijj + a_ + ^_a+), e -^0(o) = I( a+ ^_ + a_Jj + ). (2.13) 

In order to ensure that the interaction is a bosonic scalar, the a± need to be treated as anticommuting with 
the ^'s 

a±tp = —ipa±. (2-14) 

Aside from (2.14), they should be regarded as ordinary complex numbers with a± = e TMT , a + a_ = a_a + = 1. 
We can rewrite (2.11) in a slightly different form. Rearranging we have 

A f°° 

Sint = ~ / dt[(i; + e™+^_)a + S + + S-a-(^e-^ + i, + )] 
1 J-00 

A f'°° 

--± dt[(^+^ + e^)a + S+ + S-a^(ij + +ib_e- l ' T )]. (2.15) 
2 J-00 

Define the operators d and d^ 

d = S-a-, d^=a+S+, (2.16) 
which are to be treated as anticommuting with the ^>'s. The fermionized version of (2.10) becomes 

Sint = ~ / dt[{^ + e™+^_)d) + d{^e~ ia + </>+)] 

db[{il>- +^ + e ia )di + d{i)++^_e- ia )]. (2.17) 



A2 r 

2 J-< 



This is the form of the interaction we will use to calculate the reflection S-matrices. We want to calculate the 
reflection S-matrices directly from the boundary equations of motion (see below). The boundary equations 
of motion can be obtained by varying the action with respect to all the fields, including d and d) , which 
are dynamic variables, d = d(t), d) = d^(t). In order to obtain the correct boundary equations for these 
operators it is necessary to introduce the following kinetic term 

/oo 
dtd){t)d t d{t). (2.18) 
-OO 



The action which we vary is thus 



S — 5bulk + She + <Sd-frcc + Si 



'int- 



(2.19) 



The term (2.18) needs some explanation. In writing (2.18) we have considered d(t) and d^(t) to be 
one dimensional "fermionic" fields. On the other hand, (2.16) would imply that d and d) are su(2) spin 
operators. In the isotropic Kondo model it is well known that there is a screening effect [27], whereby the 
spin j operators Si act in the spin j — 1/2 representation. We claim that a similar "screening" effect also 
occurs in the anisotropic massive models. This implies that for spin 1/2 the su{2) structure is lost. The 
operators d and <$ are thus screened versions of the spin operators S±, which for spin 1/2 can be effectively 
treated as one dimensional fermions satisfying 



Equations (2.16) and (2.18) are statements of this screening. Further justification for (2.18) comes from the 
equations of motion for for S-O— and a + S + . The equations of motion can be calculated in the Heisenberg 
representation, namely d t (S-d-) = i[H, 5_a_] and d t (a+S+) = i[H,a + S + ], by making use of (2.20) and 
(2.21). This leads to the same equations of motion as from the above Lagrangian formulation in terms 
of d and d\ hence giving the same scattering matrix. Thus instead of introducing (2.16) and (2.18), we 
could have used the Hamiltonian to obtain d t (S-(i-) and d t (a + S + ), as was done in [20]. In short, (2.18) 
is consistent with a screening effect and leads to the correct equations of motion. A term similar to (2.18) 
is also introduced in the MBSG model at the free fermion point [21] and the boundary Ising model [18] to 
describe the dynamics of the boundary operators. 

Varying (2.19) with respect to all the fields, (tp±, ip±, d, d^), and requiring the boundary terms to vanish, 
we obtain six equations (at x = 0) 



d 2 =d ]2 = Q, {rf,rft} = i 



(2.20) 



and due to (2.16) anticommuting with the ip's 



dip = —ipd, d^tp = —tpd^. 



(2.21) 




(2.23b) 



(2.22b) 



(2.22a) 



(2.23a) 



idtd - -y {^+e- + V_) - y + ^+0 
id t <# = y(V>-e- iCT +?+) + y(V+ +^-e" iCT ; 



) 



(2.24b) 



(2.24a) 



7 



We see that (2.22) and (2.23) are identical, and (2.22a) and (2.22b) are adjoints, as are (2.24a) and (2.24b). 
Taking time derivatives of (2.22) and using (2.24), we get the boundary equations of motion 

ft (V+ - e-^_) = -\{X\ + A!) (V+ + e"^_) - ^AiA 2 (^_e" to +^ + ) (2.25a) 

ft (V- " e lCT ?+) = - J(A? + A^) + e^ + ) - 1axA 2 ty+e" +^_) . (2.25b) 
The bulk variation of 5buik gives the bulk equations of motion 

(ft - ft)V± + M^ ± = 0, (ft + ft)^ ± - Mip± = 0. (2.26) 

Mode expansions satisfying (2.26) take the form 

[~m~ r°° 

ip+ = J T - / d6e- l 2 {LuA-{6)e- lk - s + uA\{6y k - s ) (2.27a) 
ATF /-oo 

<P + = J— / ^e e / 2 (o;A_^)e- ife - s + w^^)e ife - s ) (2.27b) 

V_ = (V+) t , ^_ = (2.27c) 

where w = exp(i7r/4), uT = exp(— z7r/4) and k ■ x = MtcosYiQ — Mxsinh#. The variable is the rapidity and 
A±(A±) are fermion creation (annihilation) operators satisfying 

{A ± (9),Ai(6')} = 5(6-9'). (2.28) 

The boundary scattering of the particles at x — can be formulated by introducing a formal boundary 
operator B, and interpreting the boundary equations (2.25) to hold when acting on B [18]. This will give 
us a set of equations, one equation for each type of particle A\(6), of the form (sum over b) 

At(0)B = R» a (6)At(-6)B, (2.29) 

where R(9) is called the reflection S-matrix, with R b a (9) being the amplitude for an incoming particle a to 
be reflected as an outgoing particle b. If the boundary has additional structure allowing it to exist in several 
stable degenerate states, then (2.29) needs to be generalized to 

Al(9)B a = R b Ji(9)Al(-9)Bp, (2.30) 

where {B a } are the boundary states. The reflection S-matrix element R h a a {9) gives the amplitude for the 
process (a, a) — > (&,/?). Based on the "screening" discussed above, we will assume that (2.29) is sufficient 
to describe the boundary scattering involving a spin 1/2 impurity. It is known that this is sufficient in the 
massless Kondo case. However this does not entirely exclude the possible existence of boundary structure 
B a , as it may turn out that certain identifications of the matrix elements R b £ a reduce the problem to (2.29). 
That is to say, the amplitudes R b a can describe the scattering process even if there is some (degenerate) 



boundary structure. Various symmetries of the system can effectively reduce the number of independent 
amplitudes to a smaller set satisfying the same scattering constraints as R b a . Such a situation will occur 
if the soliton transition a — > b restricts the boundary transition a — > (3 to be of a specific type, with all 
amplitudes having the same soliton transition being equal. In this case the indices (a, (3) become redundant 
since the amplitudes for (a, a) — > (6, (3) are uniquely classified by the soliton quantum numbers (a, b), i.e. R h £ a 
reduces to R b a , with any boundary structure being implicitly incorporated in R b a . Poles of R h a {9) may then be 
indicative of this "hidden" boundary structure. As a specific example, suppose the boundary is described by 
two degenerate states B± with charge (+, — ) as for the solitons. Furthermore, let the symmetries (or other 
constraints) of the system be such that (i) every soliton charge (non-)conserving transition is accompanied 
by a boundary charge (non-)conserving transition and (ii) charge conjugation symmetry holds separately 
in the soliton and boundary sectors. Then R h ^ a reduces to only two independent amplitudes of the form 
R^a and R^ai where the bar denotes the charge conjugated state. The scattering process is thus effectively 
described by i?" = i?"" and J?" = R%™- (A similar reduction occurs in the bulk S-matrix at the reflectionlcss 
points.) There is also the possibility of excited "boundary bound states," arising from poles in the reflection 
amplitudes R b a {9) (see section 4). These issues will be considered in the sequel. 
For the solitons (2.29) reads 



In the case of higher spin, j > 1, it is likely that one needs to consider (2.30), with the a index taking into 
account any additional boundary structure due to the S± operators. 

We now calculate the reflection amplitudes. Substituting the mode expansions (2.27) into (2.25) gives 



Al(0)B = R^(9)Al(0)B + Rl(0)A^(0)B. 



(2.31) 





9 



(We did not explicitly write down the boundary operator B.) Solving for the reflection S-matrix according 
to (2.31) wc find 

( e p(0) + e -p(-en 

H -W- ie 2co S h(^ 7 W)V e 6 )> [2 - i6C) 

where 

e w) = coshe-^jxl + xD-^x^e 8 eP(e) = coshe-^iXj + XD-^X^e- 9 
cosh^ + 3 i 7 (A? + A|)- 5 i 7 A 1 A 2e -^ ' cosh # + 437 ( A i + A i) ~ 5a7 A * A 2 e_e 

Making use of the following relations 

e P(0) e P(-0) = e 7W e 7(-0) (2.35a) 

e 7( 4 7r/2±6) = e -7(*T/2Te) j e p(«r/2±0) = e -p(-(in/2^6)) ^ (2.35b) 

one can check that i? satisfies the unitarity and crossing relations [18] 

K(6)R b c (-0) = 6 b al ^(^-«) (2 ' 36) 

where a = —a. 

For either Ai = or X 2 = we have 

e P(0) = e 7(0) = e rt-«) = e 7(-fl) if Ai = or A 2 = 0. (2.37) 

The off-diagonal elements and J?l differ only by a phase. We will choose a such that i?7_ = Thus 
from now on we take either a = or a = n. A specific choice will be made below. Both cases can be 
obtained by applying a U(l) transformation to the fermion operators [21]. For notational convenience define 
the following 

P(9) = R+(0) = RZ(9), Q{6) = R J L{6) = R+{6). (2.38) 

The free fermion reflection amplitudes (2.33) look quite complicated. However by making a change of 
basis, their structure becomes apparent. The complex fermions can be expressed in terms of real fermions 
V>i,2 and ^1,2 

V± = -^(Vi±^2), ^± = -^=(^±#2). (2.39) 
In terms of these real fermions the action (2.19) becomes 

/oc rO />oc 

dt dx [Mdt - d x )iPi + ^(d t + d x )^ t + 2M^] - * / dt e™^^ - ^ 2 ) 
-oo J — oo 2—12 00 

/oo \ poo T poo 

dt (oiftoi + a 2 d t a 2 ) - i^- / dttyt + e^Jm-i-l dt{^ 2 -e^ 2 )a 2 , (2.40) 
-oo ^ J —oo J— oo 
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where 

Ai = A ie JCT + A 2 , A 2 = \ie la - A 2 , (2.41) 
and we have defined the operators (not to be confused with the a±) 

ai (t) = ~(^(t) d(t)), a 2 {t) = ±={d\t) + d(t)), (2.42) 



satisfying 



{ ai {t),a 2 (t)} = 0, ai (t) 2 =a 2 {tf = 1. (2.43) 



The action (2.40) is simply that of two decoupled Ising models with boundary magnetic fields. The reflection 
S-matrix for the Ising model was calculated in [18]. Note that the boundary terms are different for (^i,^) 
and (^2,^2)1 wrtn different couplings in general. The bulk and boundary equations of motion are 

(dt — d x )ipi + Mlpi — 0, (d t + d x y^ i — Mipi = (i = l,2) (2.44) 

-fttyi - e^) = ^(Vi + e^), -3 t (V>2 + e J > 2 ) = ^^2 - e*^), (at a: = 0), (2.45) 
which are just the Re and Im parts of (2.26) and (2.25). Mode expansions satisfying (2.44) are given by 

j Twr r<x> 

fo = J— / d6e- 9 / 2 (ijAi(6)e- ik - s + u>Al(6)e ik - s ) (2.46a) 

V J —00 

^. = J d6e e l 2 (uAi{6)e- ik - s + wA\{6)e ik - s ), (2.46b) 

V J — 00 

where A{ 2 (j4i j2 ) are the real fermion creation (annihilation) operators satisfying (2.28). The complex 
fermion operators can be written as 

A ± = -L(A lT iA 2 ), Al = ±J,A\±iA\). (2.47) 

Substituting (2.46) into (2.45) we can obtain the reflection S-matrices for the real fermions and hence from 
(2.47) also for the complex fermions (solitons). We find 

A{ 2 (9)B = R h2 (0)A{ 2 (-9)B (2.48) 

Ri,2{0)=Ri,2(P)f™ J) {6) (2.49) 
fli, 2 («) =itanh^Te iCT (2.50) 

A < 2 W --iisinhfl-(A li2 Te-)J' Al < 2 " 4M' (2 ' 51) 
The factor R corresponds to a fixed boundary condition for one copy of Ising fermion and a free boundary 
condition for the other copy. In the absence of any interaction, the boundary equations are (a = 0) 

^1-^=0, ^ 2 +^ 2 =0, (2.52) 
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corresponding to a free boundary for and a fixed boundary for (ip 2 ,ip 2 ). For infinite coupling 

(1^1,2 1 - * 00) the boundary conditions are 

^i+^i=0, ^ 2 -?2=0, (2.53) 

that is a fixed boundary for (^1,^1) an d a free boundary for (^2,^2)- (F° r CT = ^ (2-52) and (2.53) are 
interchanged.) So we have a flow from free to fixed (fixed to free) for ipi and ip 1 (ip 2 and tp 2 ) as the interaction 
varies from zero to infinity. In particular, the flow is different for the two Ising models. When expressed in 
terms of the complex fermions, (2.52) and (2.53) read 

V>+-VL=0, ^_-^ + = (Ai, 2 = 0) (2.54) 

V>+ + ^_ = 0, V- + = (|Ai, 2 | = 00). (2.55) 

Both equations correspond to free boundary conditions for the Dirac fermions, and apart from a phase, give 
the same reflection amplitudes. Thus there is no flow between free and fixed boundary conditions for the 
solitons. In contrast, the MBSG model has a flow between free and fixed boundary conditions as the coupling 
A goes from to 00. 

The soliton reflection amplitudes are given by 

P(6) = ±(R 1 (6) + R 2 (6)), Q(6) = ±(R 1 (6)-R 2 (6)). (2.56) 

One can check that these are the same as (2.33). In particular we have 

R 1 (6,a = ir) = R 2 (6,a = 0), R 2 (6,a = tt) = R x {e,a = 0), (2.57) 

implying, as expected from (2.33), that P{9) is independent of the choice for a, whereas Q{9) changes by an 
overall sign. 



2.1. Free Fermion Reflection S -matrix for the MK model 
We now specialize to the MK model by taking 

Ax = A, A 2 = 0. (2.58) 
Using (2.37) we get for the free fermion reflection S-matrix 

cosh(6> - 7(0)) cosh(0 — 7(0)) 

with 

~cosh0 + iA^' ~4M' [W) 
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recovering the results of [20]. From (2.56) alternative expressions are 



2 V 4 2 y Visinh^ — (A A//£ " - 1) / 2 \ 4 2 J \isinh6 - (A MK + I) 
QW = -e-itanhf4-n f ^ ^ + j*™ ~ ?> U tanh U« ■ + 



2 V 4 2/ V^in 1161 ~ ( AMK - !)/ 2 V 4 2/ V^in 1161 - ( AMX + !) 

(2.61b) 

Note that the two CDD factors in (2.61) arc different. The zero and infinite coupling limits are (the upper 
(lower) sign corresponds to A = (oo)) 

P(0) = -sech0, Q{8) = ±ie ia tanh 6, (2.62) 

which are the same (apart from phases) as the MBSG reflection amplitudes with free boundary conditions 
(A = in (1.6)). These results confirm that in the soliton basis there is no flow between free and fixed 
boundary conditions for zero and infinite coupling. This can also be seen from the boundary equations 
(2.25) with A 2 = 0. 

To check if (2.61) agrees with the Kondo model we compute the massless limit. The masslcss limit for 
right-movers is obtained by letting 8 — > 8 + a, and taking a — > oo, M — > while keeping Mq = Me a /2 held 
fixed. This gives us the massless dispersion relation E = p = M$e . Taking the massless limit we find 

R1.2 - T^ CT , A C 2 DD - - tanh (tlMl |) (2.63) 

(Q _ qMK ■ \ 
— f- ~ T j ' (2 ' 64) 

where we have defined 

M o exp(0^) = ^. (2.65) 

The massless reflection S- matrix is almost entirely due to the CDD factors /^J? - Apart from a phase, this 
is the spin 1/2 reflection S-matrix for the Kondo model [27]. The Kondo Hamiltonian (1.2) is also recovered 
in the massless limit. Thus the MK model reproduces the Kondo results at the free fermion point. 

2.2. Free Fermion Reflection S- Matrix for the mMK model 
The mMK model is obtained by setting 

Ax = V2Me*\ A 2 = V2Me^°. (2.66) 

For these values 

tsmhd — jjji^i + Ag) 

giving the reflection S-matrix 

P{8) = scch8f mMK (8), Q{6) = ie ta ta.nhdf mMK {8), (2.68) 
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where f mMK (ff) is the CDD factor 

/ mM W = " W = ( !;|nh^^M. ) > A mMK = ^(A? + Ai) = cosh20 o . (2.69) 

The CDD factor f™M«(0) i s due to the factors A C J? D (6>) in (2.49). With (2.66) wc have / CDD = ^cdd and 
jm.MK ig simply equal to -/CDD The hypcrbolic f actors i n (2.68) come from the Ri^{9) linear combinations 

^(Ri +R 2 ) = -sech(0), ^(#1 - #2) = -»e i<J tanhtf. (2.70) 

We see that the reflection amplitudes are the free MBSG amplitudes (at the free fermion point) multiplied 
by a CDD factor. The dependence on the boundary coupling is solely in the CDD factor. These elements 
are simpler than those of the MK model (2.61) where the two CDD factors are different. In contrast to the 
MK model, there is no zero coupling limit. From (2.66) it is clear that both Ai and A2 cannot vanish for 
a non-zero mass. We can take an infinite coupling limit by letting (f> — > 00 or (j) — > —00, corresponding 
respectively to (Ai, A 2 ) — > (00, 0) and (Ai, A 2 ) — > (0, 00). In either case we get 

P(6) = -sech0, Q{6) = -ie la tanhtf, (2.71) 

which are the same as (2.62). 

The massless limit is taken as before with one important difference. Because the couplings are mass 
dependent, we must also let 4>o — > 00 while keeping Ai fixed. Of course A2 goes to zero. Explicitly for the 
massless limit of f mMK we have 

f mMK [e) = iMsinh^ + Mcosh2^ _^ iM e e + M e e B MK = ( Q-B™mk ^ 

iMsinh^-M cosh 2^o iM e e - M o e B MK \ 2 4 )' 



where 0% MK is defined through 



M exp(^ MK ) - ^e 2 ^ = ^. (2.73) 



2 

The resulting reflection amplitudes are identical to (2.64) 

(Q _ nmMK ■ \ 
~2 — t) - (2J4) 

Again we recover the Kondo results. Since A 2 — > 0, the mMK Hamiltonian (1.5) goes over to the Kondo 
Hamiltonian (1.2). 



3. Reflection S-Matrices away from the Free Fermion point and Integrability 

In an integrable boundary field theory, knowledge of the particle spectrum (both the bulk spectrum and 
the boundary states) can allow one to determine the reflection S-matrix, up to CDD factors, by imposing 
the constraints of boundary factorizability (i.e. boundary Yang-Baxter (BYB) equation), unitarity and 
crossing symmetry [18], [28]. For both spin 1/2 massive Kondo models, the bulk spectrum is the sine-Gordon 
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spectrum. Also as mentioned above, we will not explicitly consider any boundary state structure as in (2.30). 
Thus if the massive Kondo models are integrable their reflection S-matrices must satisfy the same constraints 
as for the MBSG model. 

The general solution to the boundary scattering constraints for a sine-Gordon bulk spectrum, at arbitrary 
(3, was described in GZ [18]. The GZ solution depends on two formal parameters. The difference between 
the MBSG, MK and mMK reflection amplitudes will then be in the relation of these formal parameters to 
the physical parameters (mass, couplings and (3). If the results of section 2 are extendible to arbitrary (3, 
then the free fermion matrices should agree with those of GZ at (3 = V&r. 

Before making this comparison let us review GZ's solution. Define 



For arbitrary A > 0, GZ's solution to the scattering equations is 

P GZ {9) ee R+(9) = cos(<e - iK6)r(6), 

P GZ (9) ee RZ{6) = cos(£ + ik6)r{6), 
Qf-{6) ee R-(6) = Q GZ {9) ee R+(0) ee Q GZ (6) = ~ sm{2tA9)r(9). 



The function r{9) can be written as 



r(0) =r o (0)n(0). 



(3.1) 

(3.2a) 

(3.2b) 
(3.2c) 

(3.3) 



The factor r (9) is independent of the boundary couplings and ensures that the crossing symmetry constraint 
is satisfied. The boundary dependence is contained in ri(#), which takes the form 



n(9) 



cos £ 



(3.4) 



Here a{x,9) is a function not to be confused with the phase a. The variables k,t;,r) and $ are formal 
parameters with no dependence and satisfy 



cos 77 cosh = — i cos £ , cos 2 n + cosh 2 # = 1 + . 

k h 



(3.5) 



Thus only two of fc,£,r/ and d are independent. These parameters have to be related to the physical 
parameters, which is often the hardest part. The functions ro(6) and a(x,9) can be expressed as infinite 
products of T- functions (see [18]). Making use of [29] 

<-+°° dt ' 



lnr(z) 



f 

Jo 



t 



(z-l)e- 



1 



Re(z) > 0, 



the infinite T products can be shown to have integral representations, giving for r and a 

sinh (A + | x) y 



a{x,6) = 



cosh(A# + ix) 



exp 



j — ( 



dy sm(2A9y/n) 



2y cosh Ay sinh y 



(3.6) 



(3.7) 
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r o (0) = cxp 



, , sinh | Aw sinh (^r-) y 
dy sin(2A6»y/7r) 



y sinh | sinh 2ky 
We will take £ = nir and write for the diagonal amplitude 



(3.8) 



P^(9) = P^{9) = P^(6) = cos(£) cos(iA6)r(6). (3.9) 

At the free fermion point, A = 1, the expressions for the amplitudes simplify greatly 

pGZ/«\ (cos £ cosh 9)/k cos r] cosh cosh 9 ^ 

(— i sinh + cos 77) (— i sinh 8 + cosh 1?) (—i sinh + cos 77) {—i sinh + cosh d) 

Q GZ, 9 , = j sinh 9 cosh g 

(— i sinh + cos 77) (— i sinh 9 + cosh 1?) 



In particular, if 



or equivalcntly 



we have (for either sign) 



cos T) = — cosh 1? = i 1 , (3.12) 
fc = cos£ = ±l, (3.13) 



P f ? e z e (0) = secM, QfZjO) = -*tanh#, (3.14) 
which are the MBSG amplitudes for free boundary conditions. 

3.1. The MK Reflection S-Matrix at the Reflectionless points 

In this subsection we propose a natural extension of the MK reflection S-matrix away from the free 
fermion point. The MK reflection S-matrix (2.59) at the free fermion point can be rewritten as 

cosh 7 (fl) cosh(fl) / »sinh+(A M ^-l) \ 
1 ' (-isinh6» + l)Hsinh6> + cosh 7 (6>)) V* sinh ~(A MK - I) J [ ' 



= -je^smhjg)coshjg) / zsinh+(A M ^ - 1, , 

Vl J (-i sinh 6> + l)(-i sinh + cosh 7 (6>)) V* sinh -(A MK - 1) ' 
Comparing these with (3.10) and (3.11) we see that 



P(9) = P GZ (9, V = 0,0 = ^(9)) x/ CDD (0), Q(0) = - e -Q GZ (M = O,i? = 7 (0))x/ CDD (0), (3.16) 
where 

CDD _ zsinh+(A A ^-l) 
^ isinh-(A M ^-l)' 1 ' 

The phase difference between P{9) and is irrelevant, since if (P, Q) satisfy the scattering equations then 

so do (P, —Q). We will set a = tt for the MK model. Up to a CDD factor, we see that the MK reflection 

amplitudes agree with those of GZ provided we take r\ — and more importantly, allow $ to have the 9 

dependence d = 7(6*). However, such a 9 dependence is not allowed by the BYB constraint [18]. Unitarity 
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and crossing symmetry alone do not rule out the possibility of a dependence for the parameters (£, k, 77, 
If there is no BYB constraint, as is the case at the free fermion point, and more generally at the reflectionlcss 
points, then (3.15) is a possible solution to the scattering equations. Away from the reflectionless points, 
the BYB equation forces the scattering amplitudes to take the general form of GZ's solution (3.2) (up to 
CDD factors), with the formal parameters 9 independent. This means that (3.15) cannot be extended to 
arbitrary values of A in a manner consistent with (3.2). Since GZ's solution is the unique solution to the 
boundary scattering equations for a sine-Gordon spectrum at arbitrary A, we conclude that assuming there 
are no boundary degrees of freedom, the MK model is not integrable (i.e. there is no reflection S-matrix 
consistent with BYB, unitarity and crossing symmetry) for generic A. Again we would like to emphasize 
that for higher spin j > 1, one probably needs to solve a more general BYB equation with boundary spin 
states and amplitudes R b a a {0). In which case a comparison with the MBSG model cannot be made. 

The reflectionless points occur at integer A, for which the BYB constraint is trivially satisfied. The 
reason for this being that if A = n > 1, the bulk scattering matrix takes on a simple form with only one 
independent scattering amplitude. In general the bulk sine-Gordon scattering matrix is 

a(6) = S++(9) = SZZ(O) = sin[A(7r + i6)]Z(8) (3.18a) 

b(6) = S+Z(6) = SZ+{6) =-sm(iA6)Z(6) (3.18b) 
c(6) = SZ1(6) = SZZ(6) = sm(An)Z(9), (3.18c) 
where the function Z(6) can be found in [30]. For integer A we have 

6(6>) = a(6) (A odd) or b{6) = -a{6) (A even), and c{6) = 0. (3.19) 

The expression for a{9) reduces to 

A-l 

<0) = (-1) A II /»/a(«). (3- 2 °) 



n=l 

where 

sinh {\\9 + ina\) 

fa = ■ rrfe ■ ■ ( 3 - 21 ) 

sinh [^[o — iira\) 

One can check that for (3.19) there is no constraint from the BYB equation. Since the off-diagonal bulk 
amplitude c(6) vanishes, these integer points are referred to as the reflectionless points. 

With no BYB constraint, it should be possible to extend (3.15) to the reflectionless points. In terms of 
i?i,2, the unitarity and crossing relations for integer A become 

Unitarity: 

R^R^-B) = 1, R 2 {6)R 2 {-6) = 1. (3.22) 

Crossing: 
Odd A 

R^i-k/2-6) =a(26)R 1 (iir/2 + 6), R 2 (in/2 - 6) = a{26) R 2 (m/2 + 6) (3.23a) 
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Even A 

R^in/2 -&) = -a{26) R 2 {in/2 + 6), R 2 {in/2 -6) = -a{26) R^in/2 + 6). (3.23b) 



Since a(— 6) = l/a(6), the two crossing equations for even A are identical. 

Now recall the free fermion reflection amplitudes in the Ising notation (2.56), where (<r = n) 

Ri,2(e)=Ri,2(6)f^ D (e) (3.24) 
^ ,„n / 7r \ , rnD ( i sinh 9 + (A MK ± 1) \ . . 

Considering the structure of the reflection matrices (3.24) for A = 1, we will look for solutions of (3.22) and 
(3.23) of the following form 

RM = h(0)f? DD (O), R 2 {6) = r~ 2 (0)/ 2 CDD (0), (3.26) 

where f\, 2 {0) are independent of the boundary coupling A. The boundary dependence being contained in the 
CDD factors. For odd A, the CDD factors can be different for R\ and R 2 , as they are for A = 1. However 
for even A, since a(29) is independent of A, the crossing relation implies that the CDD factors must be the 
same. Requiring the CDD factors to agree with (3.25) for A = 1 suggests we take the general form 

/■cdd _ ^sinhg + r 1 , 2 (A,A 2 /G) 

Jl < 2 ~ isinh0-r li2 (A,A 2 /G)' (6 <} 

where ri j2 (A, A 2 /G) are functions of the dimensionless parameters A and A 2 /G, satisfying 

r li2 ^A=l,^ = A M *±1 (3.28) 

A 2 \ _ /. A 2 



Y l (A = 2n, - = T 2 (A = 2n, - . (3.29) 

Note that (3.28) makes sense since at the free fermion point M oc G. 

The crossing relations that need to be satisfied can now be written as 

a(26) fi, 2 (6> + in/2)f h2 {6 - in/2) = 1 (A odd) (3.30a) 

-a(20) h{0 + in/2)r 2 {6 - in/2) = 1 (A even) (3.30b) 
To solve these equations we make use of the following expression for a(20) 

A 

a(28) = H /_ n/2A (0 - in/2)f_ n/2A (6 + in/2), (3.31) 

71=1 

where the / a 's are given by (3.21). Since 

h, 2 (6)=R 1 (6), (3.32) 
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we take f\{9) to be 

A 

fi(0) = J] /"/2a(*)- (3-33) 
This will satisfy the odd A crossing relation for f\ (9) . Similarly, knowing that 

-fy 2 (0-in) = R 2 (e), (3.34) 

we take for f 2 (9) 

A 

?2(8) = -Hfn/2A(6-iiT)- (3.35) 

n=l 

Using 

f a {0 + i2irm) = f a+2m (0) = f a (6), (3.36) 

where m is any integer, we see that (3.35) satisfies the odd A crossing relation for r 2 {6). Unitarity for f\_ 2 
is also easily checked with the relations 

f a (0)f a (-e) = 1, f a {-6) = f_ a {0). (3.37) 

The solutions (3.33) and (3.35) are not valid for even A. These expressions imply 

A 

h(9 + in/2)r 2 (6 - nr/2) = - J] f 2 n/2K {6 + in/2), (3.38) 

n=l 

which does not satisfy the crossing relation (3.30b). Possible solutions for (3.30b) are 

A 

h(9) = -r 2 (9) = n/n/2A(0), (3-39) 

71=1 

and 

A 

h{9) = -h{9) = [] /n/2A(A - iTr). (3.40) 
n=l 

We can now combine the odd and even A solutions into a single expression. For any integer A we have 

A 

h(9) = 1] /™/2A (° ~ ^Pi(A)) (3.41) 

n=l 
A 

r 2 (0) = - J] (* - i*P*W) (3.42) 



n=l 

where for odd A 

pi(A = 2n+l) = 0, p 2 (A = 2n+l) = 1, (3.43) 
and for even A we can have two choices 

pi(A = 2n) = 0, p 2 (A = 2n) = 0, (3.44a) 
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or 

pi(A = 2n) = l, p 2 (A = 2n) = l. (3.44b) 

These expressions for pi^(A) arc mod 2 because of (3.36). If it were not for the specific forms (3.25) for 
A = 1, we would also have two choices for odd A. As an example, we can take for (3.44a) 

pi (A) = 0, pa(A) = sin 2 (AV2). (3.45) 

Putting everything together we have for any integer A 

D ... Asinh6> + ri, 2 (A,A 2 /G)\ -A- , . .... .„ 

RiAO) = T ( lsinhe _ ri2(AX2/G) ) II /-/»a (' - *TPi, a (A)) , (3-46) 



71=1 



giving for the soliton reflection S-matrix at the reflectionless points 

A 



«•> He f£S3$$&) n /./» « - -«(a» . 



— 2 vZsinhe-r.^AVG);^ 
Note that for even A 

P{9) = 0, A even. (3.48) 

Further restrictions on i\ 2 (A, A 2 /G) and pi j2 (A) can be obtained by analyzing the pole structure (see section 
4), and in addition I\ 2 (A, A 2 /G) is also constrained by the massless limit. 

The massless limit for arbitrary A is obtained by taking G — > 0, along with 9 — > # + a and a — > oo, while 
keeping G - G( A+1 )/ 2A e a /2 held fixed. In general M cx G( A+1 )/ 2A , thus G( A+1 >/ 2A acts as a mass, with G 
replacing M away from the free fermion point. This gives our "massless limit prescription". In order for 
the massless limit to make sense and be of the form (2.64), the limits limc^o G( A+1 )/ 2A Ti j2 (A, A 2 /G) must 
be well-defined and equal. Considering (3.28), a reasonable choice for I\ 2 is to take 

/ x2\ /A 2s(A+l)/2A 

r li2 (A,-J= c (A)^-j +r li2 (A), (3.49) 

where c(A) and n j2 (A) are functions chosen to satisfy (3.28) and (3.29). The values of n j2 (A) are irrelevant 
in the massless limit. Following this procedure we find 

— iie^-D, f c^ {e) ^_ tanh (^l^l_^ (3.50) 

(o _ qMK ■ \ 
— 2 — t)' (3,51) 

where 9g K is A dependent, defined to satisfy 

Go cxp 9™ K = Hm G( A+1 )/ 2A r! (a, = hm G ( A+1 >/ 2A r 2 (a, . (3.52) 
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We recover the Kondo reflection S-matrix for all integer values of A. Because ri^{6) is independent of the 
boundary, it only contributes a phase. This phase is actually necessary to satisfy the massless crossing 
relation. In the massless unitarity relation the phase cancels since it appears conjugated in Q{—9). Massless 
unitarity and crossing relations are discussed in [31]. 

3.2. The mMK Reflection S-Matrix for general /3 

The mMK scattering amplitudes at the free fermion point are easily related to GZ's general solution. 
Comparing (2.68) and (3.14) we find 

P(9) = Pg z Mf mMK (6), Q{6) = Qil{6)f nMK {e), (3.53) 

where we have set the phase a = it. These are just the MBSG amplitudes with the free boundary condition 
multiplied by a CDD factor, with the formal parameters satisfying (3.12) and (3.13). Without loss of 
generality, we will fix £ to be for all A. For real (3.12) gives the values 

r}(A = 1) = 7T => fc(A = 1) = 1, i?(A = l) = 0. (3.54) 

Unlike the MK model, (3.53) and (3.54) are consistent with the generic result (3.2), thus allowing for an 
extension to arbitrary A. Taking into account the structure of the free fermion amplitudes (3.53), we propose 
the following soliton reflection S-matrix for the mMK model 

P(0) = P flcc (6)f mMK (9), Q{6) = Q liee (6)f mMK (6), (3.55) 

where Pf ree (#) and Qf vec (0) are the MBSG amplitudes for free boundary conditions 



A 2 

P free (0) = cosh(A0)r o {9)a [-j^-,0) *(0, 9) (3.56a) 



QfreeW = ~% Smh(Ag) ^^L ((?)(7 (^,0) *(0,6), (3.56b) 



sin(A^/2) uw \f3 2 
and 

fmMK = ismh9 + T(A,j) ) 
ismhO -T(A,(j) y 

with r(A, (f>o) some function satisfying 

r(A = l,<£ ) = cosh2(£o. (3.58) 

The CDD factor (3.57) is a generalization of (2.69) to arbitrary values of A. In (3.56) the parameter rj has 
been given a A dependence 

^ = ^ A+1 ) = ^^ fc ( A ) = sin(AW' ^ 
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while #(A) = i9(A = 1) = 0. This serves to produce the correct pole structure, namely a simple pole at 
9 = iir/2 for all A [18]. It appears that Qf ree is singular for even A because of the sin(A7r/2) factor. But an 
identical factor appears in the numerator of a(4n 2 /P 2 , 9) (3.7), making (3.56b) well-defined for all A. Since 
at the free fermion point A = \J2M oc \/G, and G is the only available bulk parameter, we must have A 2 oc G 
for arbitrary A, thus explaining the relation (1.5c). 

The main feature of (3.55) is that the boundary dependence occurs only in the CDD factor, just as for 
MK model. The minimal (non-CDD) parts Pf ree and Qhce are independent of </>o- In contrast, the boundary 
dependence for the MBSG model is in general more complicated, with the formal parameters r\ and $ 
(and hence the factor r\(9)) being functions of the physical parameters. With the boundary dependence 
isolated, the massless limit of the minimal factors follows by taking 9 — * oo. The exact dependence of 
r] and i? on A is not important in the massless limit as long as they are independent of the boundary 
coupling. The massless limit for the CDD factor is taken in the same manner as for the MK model, with the 
difference that </> also has to be appropriately scaled such that \fGe^ a is held fixed. Holding \fGe^ a fixed 
amounts to keeping the coupling Ai constant as G — > 0. The massless limit restricts T(A,(j) ) to the extent 
that lim G ^ o < ^ o ^ oo G( A + 1 )/ 2A r(A,0o), with y/Ge't"' fixed, must be well-defined. Making use of the integral 
representations (3.7) and (3.8) we find for the functions r (9) and a(x,6) 

lim r a (9) = 1, lim a{x,9) = 2 cos(x)e" ra lim e~ Ae , (3.60) 

which gives for the massless amplitudes 

P(9) = 0, Q(9) = ie-*i ( A+1 > tanh ~ - ^ , (3.61) 

where 9 B lMK as a function of A is defined through 

G exp^ MK = lim G( A+1 )/ 2A r(A,0 o ) . . (3.62) 

G^O, 0o^oo VGe^o fixed 

So up to some phase we obtain the correct Kondo results for all A. This suggests that the integrable massive 
generalization of Kondo is given by the mMK model (1.5). We again point out that the situation is different 
here as compared with the MBSG model. For MBSG, rj and d are boundary scale dependent. To get the 
massless limit, these parameters also have to be rescaled such that g(9, 77, i?) — > 6 — 6b remains finite, where 
g is some function giving the massless limit prescription [16]. In our case the prescription (for the minimal 
part) is very simple, let 9 — > 00. 



4. Boundary Bound States 

Important information on the spectrum of a theory can be obtained from the poles of the reflection 
amplitudes. Poles in the reflection S-matrix R can be of two types. If there are bulk bound states interacting 
with the boundary, then poles will be found at 6 = in/ 2 and 9 = z7r/2 — 9b, where 9b is the bound state pole 
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in the bulk S-matrix. These bulk states will give zero-momentum one-particle contributions to the boundary 
state \B). Secondly, poles associated with boundary bound states can also occur. These are excitations of 
the boundary ground state appearing in the region < 9 < in j2. If the pole occurs at 6 = in/2, then there 
will be a degeneracy in the ground state and again a zero-momentum state will be found in \B). At the free 
fermion point we only expect boundary bound states since there are no bulk bound states (the bulk S-matrix 
is -1). 

4-1. Boundary bound states for the MK model 

At the free fermion point we can map the MK model to the two Ising models, called Isingl for (ipi,^-^) 

and Ising2 for (^2,^2)) to study the poles. These Ising models differ in their boundary conditions at zero 

and infinite coupling. The poles of the MK reflection S-matrix occur at 

7T A 2 
61 =i-, and sinu 2 = 1 - (0 2 =iv 2 ), (4.1) 

where 9\ is due to Isingl and 9 2 = iv 2 , which only exists in the physical region if A 2 < 4M, is due to Ising2. 
As in [18] these have the following interpretation. Each Ising model has two degenerate boundary ground 
states |0,±)^ {i — 1,2}, giving two different expectation values for the spin field a(x), (<j(x)) ± cx ±1. For 
the Ising2 piece, the degeneracy is broken by the boundary magnetic field 

\0,±f B ^\0) 2 B ,\l) 2 B , (4.2) 

and the pole at #2 corresponds to the boundary bound state \1} B with energy 

E{ 2) -£7< 2) =Mcosv 2 . (4.3) 

In the Isingl case the degeneracy persists for a non-zero boundary interaction (A 7^ 0). Both states have the 
same energy 

eW=eW. (4.4) 

Thus for A 2 < AM we have two sets of degenerate states 

{|0,+>>> 2 S ,|0,-)>> 2 B } and {|0, +^11)^10,-^11)^}, (4.5) 

associated with the poles at 9\ and 9 2 respectively. This suggests that one should actually regard the 
boundary as having structure B a , where a = ± corresponds to the above degenerate states. Because of the 
pole at in/2, there will be a zero- momentum contribution to \B) of the form A\(0) = -^(A+(0) + A^_(0). If 
A 2 > AM, the pole at 9 2 leaves the physical strip and we are left with the degenerate ground states. 

It is not clear what this pole structure implies for the states of 4> and the spin operators S± . We make 
the following remark. For a spin 1/2 impurity along with the screening effect, the boundary is described by 
the Hilbert space 

I®i = oei. (4.6) 
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For an infinite interaction (A = oo) the boundary would exist in the | j, m) = |0,0) state. With a finite 
interaction (A 2 < 4M), the space (4.6) decomposes into two sets of degenerate states, which are some linear 
combinations of |0,0), |1,0) and |1,±1). These have to be related to the Ising states. 

Away from the free fermion point the bulk S-matrix has soliton bound states. These are the charge 
neutral breathers with masses 

/ ft 7T \ 

m„ = 2Msin(^ — J , n=l,2,...<A. (4.7) 

The bound states should give rise to poles in the reflection amplitudes. For odd A > 1, the poles in P{9) 
(3.47a) at 

TL7T 

6 = i 2A' n=1 < 2 '-- <A ' ( 48 ) 
are associated with the breathers. These poles also occur in Q(9), implying that the boundary should carry 
a (+, — ) "topological charge" since the breathers are charge neutral. This would agree with a description 
of the boundary in terms of B±, with the degenerate states (4.5) having different charges. To account for 
the breathers at even A > 1 we must take (3.44a) for Pi(A), which will reproduce the poles (4.8). The 
diagonal amplitude P(9) vanishes for even A. This would make sense if at even A the boundary charge is 
not conserved, leading to a change in the soliton charge. 

With the choice (3.44a), the pole at 9i = i-x/2 is present for all integer A and the ground state is 
twofold degenerate. The existence of a boundary bound state for A > 1 depends on ri !2 (A, A 2 /G). If we 
have — 1 < ri i2 (A, A 2 /G) < 0, then a pole will exist in the physical strip at 

9 2 = iv 2 , sinv 2 = -r li2 (A, A 2 /G), (4.9) 

leading to a bound state. However for A > 1, a map to the Ising picture does not exist and the boundary 
states cannot be explained as in (4.5). For odd A we can have —1 < ri j2 (A = 2n + l, A 2 /G) < and Ti ^ r 2 , 
which will lead to two distinct bound states. 

4-2. Boundary bound states for the mMK model 

The pole structure for the mMK model is identical to that of the MBSG model with free boundary 
conditions, apart from any poles due to the CDD factor. We will only comment on the boundary bound 
states. At the free fermion point the mMK model can be written as a sum of two Ising models, Isingl and 
Ising2, with couplings Ai and A 2 . Consider for the moment the more general system (2.40). The first Ising 
copy leads to the pole at 9\ — in/ 2 as for the MK model. The second copy contains the CDD factor 

/•cdd = _ / »sinhfl+(A 2 -l) \ 

h ~ \zsinh# - (A 2 - 1)/ ' ( j 

which gives a boundary bound state provided the "magnetic field" is not too strong 



A 2 < 4M. (4.11) 
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In terms of Ai and A 2 (<r = ir) 

A| = A? + A^ + 2A 2 A 2 . (4.12) 

Now specialize to the mMK system with 

Ai, 2 = \/2Me ± ^°, (4.13) 

giving 

\ 2 2 = 4M(cosh20 o + 1) > 4M. (4.14) 

Thus there is no boundary bound state because the coupling | A2 1 is too large for all values of <j>Q. The poles 
for the CDD factor f™ MK {6) 

f mMK {e) = _ /a GDD ((?) = ( * sinh g + COsh 20 o \ ^ 

\ i sinh 8 — cosh 20 o / 



occur outside the physical strip at 



fl = ±2^o-i|. (4.16) 



So the boundary states at A = 1 consist of two ground states 

{|0,+»)|,|0,-)»|}, (4.17) 

corresponding to the pole at z7r/2. 

By choosing 77 = |(A + 1), we maintain the pole at in/ 2 for all A, though its interpretation in terms of 
the Ising states is lost for A > 1. Away from the free fermion point f mMK {Q^ has a pole at 

sinh 2 = -ir(A,(£o). (4.18) 

If —1 < T(A, 0o ) < 0, this pole will lie in the physical strip implying a bound state. 



5. Conclusions 



We have presented two massive versions of the anisotropic spin 1/2 Kondo model and discussed their 
intcgrability. Both models correspond to a sine-Gordon theory in the bulk, but differ slightly in their 
boundary interactions. Our arguments for intcgrability were based on the assumption that there arc no 
explicit boundary degrees of freedom in the scattering description, which is supported by the screening at 
the free fermion point, on constraints from the boundary Yang-Baxter equation, and comparison with the 
massive boundary sine-Gordon model at the free fermion point. We argued that the most natural massive 
extension of the Kondo model (MK) is not integrable except perhaps at the reflectionless points, and we 
proposed a reflection S-matrix for these points. The modified massive Kondo model (mMK), on the other 
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hand, we argued is integrable for all (3 and proposed an exact reflection S-matrix, which is a CDD factor 
times the massive boundary sine-Gordon solution corresponding to a free boundary condition. 

Since we did not prove the integrability of our models by studying integrals of motion, our proposed 
S-matrices and the claim that the modified massive Kondo model is integrable remain conjectures. However, 
the fact that we were able to construct an exact S-matrix, starting with an explicit calculation at the free 
fermion point and consistently (i.e. agreeing with the scattering constraints) extending the result to arbitrary 
(3, supports the integrability claim. A more concrete analysis would involve trying to construct quantum 
conserved charges for arbitrary (3. Compared to the massive boundary sine-Gordon model and the (massless) 
Kondo model, such an analysis for the massive Kondo models is more complicated due to the presence of both 
(i) interactions in the bulk and at the boundary and (ii) spin operators S± at the boundary. In relation to 
this, it would be interesting to develop conformal perturbation theory to study integrability in the presence 
of boundary operators such as S± . Lastly, a classical study of the massive models can also be attempted. 
Classical integrability can be investigated through a modified Lax pair approach as introduced in [32] to 
study boundary Toda theories. We plan to address these questions in the future. 
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Appendix 



We briefly discuss how the anisotropic spin 1/2 Kondo model can be mapped onto the one dimensional 
bosonic field theory 



HK = \f dx {(Mf + (d^f) + A (s +e ^°>/ 2 + S_e-*. 



/30(O)/2\ 



(Al) 



where A and are free parameters explained below. A more complete derivation can be found in [33] (see 
also [34]), which we will follow closely. The map involves the technique of abelian bosonization, as was 
first applied to the Kondo model by Schotte [35]. For further information on bosonization see [36]- [39] and 
references therein. 

We begin with the anisotropic spin 1/2 Kondo Hamiltonian 

H K = ^ + H int (A2) 

ka 

ffmt = E 4>*)-' c ^< + j ~y E + vi(°*)^<w) • ( A4 ) 

kk'aa-' kk'aa' 

The interaction is taken to be pointlike with only s-wave scattering. Thus the plane wave states can be 
expanded in spherical waves about the impurity, with only the angular momentum I = m = modes being 
retained. Also since we are usually interested in low energy excitations near the fermi surface, e(k) can be 
expanded as 

e(k) = e F + v F (k - k F ), (A5) 

where k — \k\. With these approximations, and measuring the energy and momentum relative to the fermi 
surface, the free part becomes (with v F = 1) 

H o = E pc ^ c p ct ' ( A6 ) 

pa 

where c^ CT create s-wave electrons of momentum k = k F +p. As is common for bosonization, in fact necessary 
[36] , we allow p to be unbounded in Hq and regularize by introducing an exponential cutoff in the interaction. 
The ground state consists of all states with p < filled. Using the fcrmion fields 

where L is the normalization length and p = 2im/L, n = 0, ±1, ±2, . . ., combined with a redefinition of the 
coupling constants, the interaction takes the form 

Hint = J ~fS z (^ T (0) - ^(0)) + ^ (s+V^tW + 5-^(0)) , (A8) 
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where S z — o z and S± = \(p x ± ia v ). The coupling constants are now dimensionless. 

To bosonize, it is convenient to separate the so-called fermion charge and spin degrees of freedom. Define 
the following charge and spin density operators for k > 

Pc(k) = ^( c P+fcT c Pf +c I+H c p^> Pc(-k)=pt(k) (A9) 
p 

Ps(k) = ^2(cl +k1 c P i-cl +kl c pl ), Ps{-k) = p\{k). (A10) 



Then the operators 



4 = J j^(-*0, a' k = J j^Ps(-k), (All) 



obey canonical boson commutation relations 

[a c k ,a$] = [a' k ,a$]=5 kk ,, (A12) 
when acting on the ground state. Furthermore, instead of Hq we can use the boson Hamiltonian 

fe>0 

to obtain the same dynamics. The charge and spin sectors separate in the linear approximation. The J z 
term can now be expressed as 

^(4M0)-4M0)) = : jS z J2e- ak/2 J^ z (a s k + a^). (A14) 

fc>0 

Here we have introduced a cutoff which eliminates values of k > a -1 ~ kp, thus regularizing divergent 
momentum sums. To write the J± interaction in terms of the boson degrees of freedom, we make use of the 
"bosonization" rule 

i>1,l(x) = -jLe-^'iW, (A15) 
V27ra 



where the fields 0fj arc 

'2tt 

fc>0 

with the 6 operators related to the charge and spin operators 



h,i = i Yl e ~ ak/2 Viz ( b ^-^ lkx b li,i elkx ) ' #J = *U (A16) 



6fcT = ^(Ofc + o*). &fcl = ^K-o*)- ( A17 ) 

For the same spin index a (T,|), the exponentials in (A15) behave as fermion fields. The only problem 
with exp(±i(j) a ) is that for different spins a ^ a' these commute. This can be fixed by using more complex 
representations, such as including Klein factors [36] (see also [39], [24]). The Klein factors serve to produce 
the correct commutation relations as well as create and destroy electrons, just as the fields and tp. The 
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exponentials alone cannot change the number of electrons. We will not consider these important issues here. 
With (A15) the J± term can be written as 

^ (S + ^ T (0) + S_$^(0)) = ^ (s+e'^W + S_e-^°)) , (A18) 

where 

= 4^ (M*) - M*j) = -iJ2 e ~ ak/2 ^u( a k e ~ ikx - a l^ lkx ) (A19) 



fe>0 

We see that the exchange interaction involves only the spin sector, the charge sector is completely decoupled. 
We will only be interested in the spin sector and drop the s superscript. Combining everything we have 

H = V ka{a k ~ ±*S x d x m + ^- (s+e^*™ + S_ e "^(o)) . (A20) 

Now we can write down the boson field theory Hamiltonian associated with (A20) . Consider the following 
action on the half line 



S=\f dx ((d t <f>f - {d x 4>f) . (A21) 

-co 

Varying the action and requiring the boundary terms to vanish gives the boundary condition 

d x <j> = 0, at x = 0. (A22) 
A mode expansion satisfying the bulk equation of motion, (d t 2 — d x 2 )(f> = 0, can be written as 

-co^^f^"^- ^)' (A23) 

where kx = \k\t — kx. The boundary condition (A22) requires that a(k) = a(—k). This allows us to separate 
(f> into its left and right moving parts as follows 

4> = H + 4>r, HAt, x) = ~i r "S= -7kt Uk)e- lk(t±x) - a\k)e^ t±x A . (A24) 
Since <Pl{x) — <Pr{— x), the free Hamiltonian obtained from (A21) is 

/>oo 

= / dkka*(k)a(k), (A25) 
Jo 

which is just the continuum limit of the free part of (A20). Also taking the continuum limit of <j> we find 
that 

0(t, x = 0) = 0(t, 0) = 2<f> L (t, 0) = 2fo(*, 0), dj(t, x = 0) = d t <j>(t, 0) = n(t, 0). (A26) 

Of course we should include the exponential cutoff in (A23). This shows that the anisotropic Kondo model 
(to be exact, the spin sector) is equivalent to the scalar boundary field theory 

H K = \f dx {(d^f + (dAf) £_Ls z 7r(0) + J~~ (s + e^^ + S _e- i ^°A . (A27) 
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Finally to obtain (Al) we make the unitary transformation 



where /i is a constant. The transformed Hamiltonian, apart from constants, is found to be 



2j_oo V 4 v27r / 47ra V / 

(A29) 

Choosing /j, to be 

"— T^r- < A30) 

the J z term cancels, leaving the boundary field theory Hamiltonian (Al) 

H K = \f dx ((Mf + (d x <P) 2 ) + A (s+e^' 2 + S_ e - 4 W°>/ 2 ) , (A31) 

2j-oo ^ ' 



A = -^, /3 = V27f2-^). (A32) 



where 

A = . _ 

47ra \ 27r / 

The coupling A has a canonical dimension of mass whereas [3 is dimensionless. For an antiferromagnetic J Zl 
47r > J z > 0, we have < (3 < V87T- If we take J± to be small, then the isotropic Kondo model corresponds 
to J z being small J z ~ Jj. . Thus the isotropic point occurs approximately at /3 = V87T. We have shown that 
the anisotropic Kondo model can be described by the bosonic boundary field theory (A31). 
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